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Abstract 

Yukawa potentials may be long ranged when the Debye screening length is large. In 
computer simulations, such long ranged potentials have to be taken into account with 
convenient algorithms to avoid systematic bias in the sampling of the phase space. Re- 
cently, we have provided Ewald sums for quasi-two dimensional systems with Yukawa 
interaction potentials [M. Mazars, J. Chem. Phys., 126, 056101 (2007) and M. Mazars, 
Mol. Phys., Paper 1]. Sometimes, Lekner sums are used as an alternative to Ewald sums 
for Coulomb systems. In the present work, we derive the Lekner sums for quasi-two di- 
mensional systems with Yukawa interaction potentials and we give some numerical tests 
for pratical implementations. The main result of this paper is to outline that Lekner 
sums cannot be considered as an alternative to Ewald sums for Yukawa potentials. As 
a conclusion to this work : Lekner sums should not be used for quasi-two dimensional 
systems with Yukawa interaction potentials. 
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1 Introduction 



This work follows a previous paper on the derivation of Ewald sums for Yukawa potentials 
in quasi-two-dimensional systems [U |2] . The concern of the present work is to derive the 
Lekner sums for Yukawa potentials and to explore their possible applicabilities to quasi- 
two-dimensional systems. 

J. Lekner has developed his method as an alternative to Ewald sums for Coulomb poten- 
tial in systems with periodic boundary conditions [3], H] ; from an analytical point of view 
Lekner and Ewald sums are fully equivalent. In Lekner sums, the summations over the 
periodic images are transformed as summations over modified Bessel functions Kq (or Ki 
for molecular dynamics implementations) ; because of the asymptotic behaviour for large 
values of the argument of Bessel functions, the summations are rapidly convergent ; but, 
because of the singular behaviour of Bessel functions for small values of the argument, the 
summations converge very slowly for some particular configurations of the pair of particles 
[HI [6]. This last property of Lekner sums introduces some complications in numerical im- 
plementations of Lekner sums. There are mainly two methods to improve the convergence 
rate of these summations : the Lekner-cyclic [4], [7] and Lekner-Sperb [HI El [ID] methods. 
A review of the uses of Lekner sums for Coulomb potential in computer simulations has 
already been done in refs. O [7] and would not be reproduced here. 

In the present work, we derive Lekner sums for quasi-two dimensional systems with 
Yukawa interaction potentials and we give some numerical tests for pratical implementa- 
tions that use the Lekner-cyclic method. The paper is organised as follow. In section 2, 
we derive the Lekner sums for Yukawa potentials in quasi-two dimensional systems and 
we show from an analytical point of view that Lekner and Ewald sums are equivalent. 
In section 3, we give some numerical tests for implementations of Lekner sums, these 
tests show that the use of Lekner sums for Yukawa potential should be avoided since to 
reproduce with accuracy the Coulomb limit of Yukawa potentials, one has to modify the 
summations by using Taylor expansions for small screening parameter /t that would not 
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allow to implement Lekner sums for any value of k. The conclusion in section 4 states 
that one should use Ewald sums for Yukawa potentials in quasi-two dimensional systems 

mm- 



2 Lekner sums for Yukawa interactions in quasi-two 
dimensional systems. 

We consider a system made of particles interacting via Yukawa potentials. The simu- 
lation box have partial periodic boundary conditions, with spatial periodicity and Ly, 
applied respectively to directions x and y, whereas no periodic boundary conditions are 
taken in the third direction z, parallel to the unitary vector e^. In the simulation box, 
the position of the particle i is defined by = (xj, i/i, Zi). 
The particle-particle interaction energy is given by 

1 Af TV 1 TV 

i?,,(Yukawa; «:) = - E E QiQi^i^ij) + ^ E ^^^^o (1) 

^ i=i i^j ^ i=\ 
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^^ exp(-.|r + n|) ^ ^ .M'. I n |) 

n \ r + n \ ^^.^^ | n \ 



where we used the condensed notations 



r + n 1= ^{x + n^L^y + {y + UyLyY + z"^ and | n |= ^JnlLl + nlLl 



with and Uy integer numbers associated with periodic images of the box. 

As for the Ewald sums [1], the lattice sum is transformed by using the integral relation 
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we obtain 



X 



[ f exp(-(f +mfLjt)] 



(4) 



Then, the Poisson-Jacobi identity 

eynf — fn, -I- 

L. 



E exp(-(n + f )^L^t) = f v^[l + 2Ecos(27rpf )e-'^V/L^*] (5) 



n=— oo 

is apphed to only one of the summations over the periodic images, it gives 



thus, taking into account the integral 



jcxp[^ — -BH]= 2 Ko{AB) 



(6) 
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we find 



2 r / 1 

$(t-) = y, ^; L^, Lj,; 00, 00) = — ^ Kq [k-^/ (y + mLy)^ + 
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+ Agcos(27rpf ) f KoL (K2 + ^)((y + mL,)2 + z2)] 

p=i -t^a; m=-oo \ 
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where the notation y, 2;; L^;, Lj^; 00, 00) indicates that the summations over p and m 
are not truncated. 

The self contribution $0 is needed to compute the particle-particle energy. This contri- 
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bution reads as 
$0 = ^o{L^,Ly; 00,00) = 



exp(— K I n I) 



n 



o +00 +00 

f EE Kc 

^ p=l m=l 
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— ^ Ko(m/tLy) - — ln(l - exp(-KL^)) 

m=l 

(9) 



The particle-particle interaction energy i?cc (Yukawa; /t) is given by Eqs.(l), (8) and (9), 
for any value of the inverse screening length k. 

In the Coulomb limit (k 0), the singular terms are the first contribution of the right 
handed side of Eq.(8) and the last two contributions of Eq.(9). For the singular contribu- 
tion in Eq.(8), we find [10] 
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with A = L^Ly, and for the singular terms in Eq.(9), one has 
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ln(l - exp(-fi;L^.)) = lia{nL^) - ]-kL^ + ^{kL^Y - \{KLxf + o{k'^ 
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with 7 ^ 0.5772156649... the Euler's constant. The first expansion in (11) have been 
obtained with the help of Schlomilch series [TI] . 

On Eq.(ll), one should note that the singular terms as In/t are cancelled by adding them 
in $0- For numerical applications, Eqs.(lO) and (11) arise some difficulties in the Coulomb 
limit : expansions in the right hand side of these equations are exact as /t — 0, as long 
as the Bessel series are not truncated. However, in applications, one has to truncate 
these series to an index m-c and as k becomes smaller, rric has to be taken greater so that 



to obtain a reasonable accuracy for expansion in Eqs.(lO) and (11). This point will be 
illustrated in the next section. 
Therefore, in the Coulomb limit, we have 
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and 



*o = S + ^ i; i; K„[2.pm^5.] + ^[t + In (-,1^)] + . + o(^) (13) 
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Thus, as K ^ 0, we find 



^x 



£^cc(Yukawa; k ^ 0) = £Jcc(Coulomb) + ^( X] — ^ 
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(12) 



(14) 



with £'cc(Coulomb) the particle-particle interaction energy of a quasi-two-dimensional 
system of particles interacting with Coulomb potentials and computed with the Lekner 
sums as 
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where the prime in summations over the particles indicates that m = is not included if 
i = j. Equation (14) is exactly the same as the one found for Ewald sums, this outlines 
the consistency between the Lekner and Ewald methods for Yukawa potentials. 
In numerical uses, summations over p and m in Eqs.(8) and (9) have to be truncated to 



some finite values pc and mc, therefore, in numerical implementations of Lekner sums, 
for both Coulomb and Yukawa potentials, one should take care of the slow convergence 
of the sums over the Bessel functions for some particular configurations of the system of 
particles [3 El E] . There are mainly two ways to counter the slow convergence behaviour 
of the Lekner sums, Lekner-Sperb [10] and the Lekner-cyclic methods [7]. 
The method proposed by Sperb begins with an analytical identity similar to Eq.(8), then 
an alternative expression which converges faster than the sums over the Bessel function, 
for the configurations leading to a slow convergence behaviour is derived and used as an 
alternative to the slow convergence |H1 E] ; with this method, some cautions are needed 
to take properly into account the electroneutrality of the system [TU] . 
In the Lekner-cyclic method, one uses a symmetry property of the analytical forms of 
Lekner sums. Namely, if in Eqs.(8) and (9), one substitutes x to y and to Ly then the 
energy E'cc (Yukawa; k) has to remain unchanged (for a detailed discussion and numeri- 
cal tests about Coulomb potential, see refs.[7lE]) ; in the following, when substitutions 
X ^ y and y x are done, the new functions obtained from Eq.(8) and (9) are noted 
'i'{x,y,z;L^,Ly;pc,mc) and ^o(^x, Pc, "^c), i-e. 

^(x, y, z; L^, Ly-, pc, rric) = $(?/, x, z; Ly, L^; pc, rric) 

(16) 

"ifoiL^, Ly] Pc, rric) = ^o{Ly, L^; pc, rric) 
and obviously, we have 

lim Mx,y,Z]L^,Ly]pc,mc) = lim J/{x,y, z; L^, Ly;pc,mc) 

(Pc, rric )^ (CO, oo) (Pc,nic)^(oo,oo) 



(17) 



lim ^o{L^,Ly;pc,mc) =^ lim "^oiLx, Ly;pc,mc 

^ (pc,»Tic)^(oo,oo) (pc,nic)^(oo,oo) 



Then, to implement the Lekner-cyclic method, a criterion has to be chosen to determinate 
which function, $ or \1', has to be used to compute _Ecc (Yukawa; k) (or i^cc (Coulomb)). 
Some numerical tests to help to choose criterions for different systems will be provided 
in the next section. In the remaining of the present section we will focus on the singular 



term as k — > 0. 

In the Coulomb limit, the singular term in Eq.(14) is cancelled if the system of particles 
fulfills the electroneutrality (i.e. J2f=i Qi = 0), as in the computations done in the previous 
paper to obtain the Ewald sums [1]. Such electroneutrality implies that there is some 
attractive interactions between particles (for instance systems of macroions with opposite 
charge), but for most of systems in interest, as colloids or dust plasmas, interactions 
between the macroions are always repulsive. Thus, a neutralizing background has to 
be included in the systems to permit to obtain the correct Coulomb limit. In ref.[T], 
several kinds of neutralizing background have been considered. In the present work we 
will restrict ourselves to the monolayer neutralizing background (system (a) of ref.[T]), 
since the computation for other backgrounds may be easily done following the methods 
below and in the previous paper. 

The monolayer neutralizing background is a plan with an uniform surface charge density 
localised at z = 0, particles may be localised on both sides of the plan or their location 
may be restricted also to only one half-space. As in ref.p], the system considered is made 
by the monolayer and the particles ; for this system, the charge density in the right hand 
side of the Helmholtz equation is given by 

p{r) = Y.QAr-r.d + cx6{z) (18) 

i 

assuming that for all particles Qi = Q, the electroneutrality for this system reads as 

NQ + Aa = (19) 

Taking into account this background, the energy of the system is given by 

£'(Yukawa; k) = £'cc(Yukawa; k) + Ecs{k) + Ess{i^) (20) 

where Ecs{f^) is the interaction energy of the particles with the neutralizing background, 
while Ess{k) is the interaction energy of the neutralizing background with itself ; these 
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interaction energies are given by 
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then, with Eq.(3) and the Poisson-Jacobi identity Eq.(5), and after having performed the 
integral over the monolayer, we find 



Ecsi^^) 



Essi^) 



-2-71 — : 2^exp(-K I Zi I) 



A K 



1=1 



(22) 



TT- 



A K 

Therefore, in the small screening limit, we have 

^(Yukawa; k ^ 0) = Ecc(Coulomb) + 2tt—^ J2 \ \ (23) 

with £'cc(Coulomb) given by Eq.(15). Taking into account that the energy of Coulomb 
systems computed with Ewald sums or with Lekner sums agree very well [7j, then the 
result obtained in Eq.(23) is fully equivalent to the result obtained with Ewald sums. 
Thus, from an analytical point of view, Ewald and Lekner sums for quasi-two dimensional 
systems are fully equivalent. 



3 Numerical tests for Lekner sums. 

As outlined before for Coulomb potentials [SI El El, the uses of Lekner sums to obtain 
interaction energies (or forces in molecular dynamics implementations) have to be done 
cautiously. For some configurations of particles, the convergence of sums may be very slow 
[7] and many contributions have to be included to obtain rather accurate results and to 
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avoid bias in the sampling of the phase space [5J. For non-screened Coulomb interaction, 
the slow convergence rate stems from the first contribution in the right handed side of 
Eq.(15) ; the origin of this slow convergence rate is in the behaviour of the modified 
Bessel function Kq, or Ki in molecular dynamics computations, as their argument tends 
to zero. In Lekner sums for Yukawa potentials, the same kind of difficulties is encountered 
since a similar contribution is included in the energy (the second contribution in Eq.(8)). 
There are mainly two methods to improve the convergence rate of these summations : 
the Lekner-cyclic |H [7] and Lekner-Sperb [8[ |9l [10] methods. The method derived by 
Sperb consists in deriving an alternative formula that converges faster than the original 
expression involving the Bessel functions, while the Lekner-cyclic method consists in using 
the symmetry of Eq.(16). 

The Coulomb-like slow convergence in Lekner sums for Yukawa potentials is outlined 
below similarly to the study done in ref. [5]. We consider a pair of particles that carry 
charges Qi = +1 and (^2 = +1, in a box with dimension = Ly = 20 a where a = 1 
is a typical length scale of the system as a hard sphere diameter ; in these reduced units 
the energy at contact Ei2{r = cr) = 1. Particle 1 is located at (0,0,0) and particle 2 at 
^"12 = (3^1252/125^12)- The difference A_Ei2(/t; Pc5 "^c) is the difference between the particle- 
particle interaction energy computed by using $ or for the potential, it is defined by 

AEi2{K;pc,mc) = QiQ2{^ir 12; L^, Ly;pc,mc) - "^{vu; L^^, Ly;pc,mc)) (24) 

On Figure 1, we show a topographical representation of AEi2{k,;Pc,^) for several values 
of K and pc, such as —5 < xu < 5 and —5 < yu < 5 and zu = 0.1, while on Figure 2, we 
show a surface representation of the same quantity for same configurations of the pair. 
The simulations box dimensions are = Ly = 20 a, on Figures 1 and 2, we represent 
only situations where | x | < 5 and | ?/ | < 5 to avoid irrelevant complication induced by the 
minimum image convention that should be used to compute the energy. The oscillatory 
behaviour of AEi2{k,; Pc, rric) has the same interpretation as for the Coulomb potential [5]. 
Further numerical examinations of the analytical derivations done in the previous section- 
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show that the slow convergence rate for these contributions to Lekner sums for Yukawa 
potentials may be handled by exactly the same means as it is done for Coulomb potentials. 
Thus, to implement the Lekner-cyclic method for systems with Yukawa interaction poten- 
tials one may proceed exactly as it is done for Coulomb potentials. Similarly to Eq.(28) 
in ref.jTj, the criterion to choose which one formula between \l/ and $ (cf. Eq.(16)) has to 
be used to compute the energy of the pair, may be taken as 

^ii(Pc,"^c) for > 



Lx Ly 

$,,(p„m,) for M<M 

and, similarly to Eq.(31) of ref. [7], the truncation index pc for the sums over the cosine 
functions can be chosen such as 



M i^^ + -^PciVij^ Zij)] [{vij + mLyf + zl] > 19. 

As shown on Figures 1 and 2, the value of pc may be reduced significantly if k, is large ; 
however, as it is shown in Table 2 and Eq.(47) of ref. [Ij, for these large values of k a direct 
truncation of the potential and the minimum image convention are sufficient to obtain a 
good accuracy. Therefore, to implement the Lekner-cyclic method for Yukawa potentials, 
one may handle the Coulomb-like slow convergence behaviour of the sums exactly as it is 
done for non-screened Coulomb potentials. 

For Yukawa potentials, the Lekner sums derived in the previous section have an additional 
slow convergence behaviour : in Eq.(8), the first contribution in the right handed side 
involves a single summation over Bessel functions ; as k ^ and k 7^ 0, because of the 
asymptotic behaviour of Ko(x) as x 0, a lot of contributions in the summation have 
to be included. As we will see below, this additional slow convergence behaviour is very 
penalising for implementations of Lekner sums for Yukawa potentials. 
For a system of two Yukawa particles carrying charge Qi = +1 and Q2 = —1, the 
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Coulomb interaction energy can be obtained from Eq.(14) as k 0. Since for this 
system, we have Qi + Q2 = 0, no background is needed to restore the electroneutrahty of 
the system. Such test configurations have already been considered for Coulomb potentials 
in refs.p, [31 HJ [5], [121 [E] that provide some reference points for the Coulomb limit, for 
these configurations we take L^. = Lj^ = 1 as in original works. On Figure 3, we give 
E'cc (Yukawa; k) as functions of n computed with Eqs.(8) and (9) inserted into Eq.(l), 
for configurations already considered in numerical tests of Lekner sums for Coulomb 
potentials and in numerical tests of Ewald sums for Yukawa potential [Ij. On Fig. 3, 
a thick horizontal line, for each reference configurations, indicates the limit obtained 
with a straightforward use of Lekner sums with = 50 and rric = 4. In Table 1, we 
report, for each of five configurations already considered for the Coulomb potentials in 
refs.[3l [U [T2I [131 [S]; the Coulomb limit (k = 0) and the values obtained for Yukawa 
potentials with : k = 10~^, Pc = 50, nic = 3 and 50 [Note : In reduced units and with 
Lx = Ly = 20, a value as small as k = 10"^ for a box of side L = 1 corresponds to a 
reduced screening parameter k* = Ka = 5 x 10~^ , cf. ref.jl]]. These data show that even 
by choosing a value as large as = 50 the Coulomb limit is not well reproduced. Values 
as large as rric = 50 for the truncation of the summations over the Bessel functions are 
not a reasonable choice, moreover such choice do not allow to reproduce with a sufficient 
accuracy the Coulomb limit. 

This lost of accuracy of the Lekner sums for small k, stems from the first contribution 
in the right hand side of Eq.(8) ; it is different from the slow convergence observed in 
Lekner sums for Coulomb potentials and increasing pc, as it should be done for Coulomb 
potential, does not allow to achieve a better accuracy. Figure 4 shows that the increase 
of Pc do not permit to obtain the Coulomb limit (/t = 0.0001) with a good accuracy 
; for the configuration (xi2, ?/i2, -212) = (0.4,0.4,0.1), considered in Fig.4, any value of 
Pc > 10 give the same results for a nic fixed, but a strong dependence on tjIc, at pc fixed, 
is observed. The same is true for the other test configurations. Examination of Eq.(8) 
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allow to understand easily this independence on pc : the first contribution is independent 
of Pc- 

To overcome this lost of accuracy for small k in Lekner sums by using Eqs.(8) and (9), 
we replace the first summation in Eq.(8) by its Taylor's expansion 
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Thus, the particle-particle interaction energy is corrected as 
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More precisely, the particle-particle interaction energy, computed with the modified Lekner 
sums of Eqs.(25) and (26), is 



£;cc(Yukawa; k) = ^ E 'QiQi E E ^'^V^) Ko[ 
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This partial expansion of the particle-particle interaction energy allows to obtain the 
correct Coulomb limit as k is small, while keeping rric to a reasonable value. This is 
illustrated on Table 1 and on Figure 5, where £^cc(Yukawa; k), computed with rric — 3 and 
the modified Lekner sums, is plotted for the configuration (0.4, 0.4, 0.1). As shown on the 
inset of Figure 5, the corrected value for the energy may be obtained with accuracy for 

K < 0.1. 

For values of the screening parameter k between 0.1 and 2, that are interesting values of the 
screening parameter for many systems, the accuracy of the corrected value £'cc(Yukawa; k) 
has to be improved by adding high order contributions to and $o to keep a reasonable 
value for rric- Contributions in and are given by 

3 

RplT,ioCO^I — VtJj.UO.\\ - -- 

87r2 



C2{ry 



^ Re[Li3(n)] - ^ Re[Li.(n)] - 



6 LI 
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with Re(2;) the real part of the complex number the polylogarithm functions lAi,[z) and 
Q defined by 



Q = exp(-27r(^ + iy)lLy) and Li^(Q) = ^ —9P 

p=i P 



One has also to add to $o the contributions 



TT^ LI 
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where ({z) is the Riemann's Zeta function defined by 



OO -1 

p=l p 
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For K > 2, the particle-particle interaction energies may be evaluated with a good accuracy 
by using Eqs.(8) and (9) with = 3. From data on Figure 5 and by using the original 
Lekner sums in Eq.(8) and (9) at k = 2 and with pc = 50, we have for rric = 3, Ecc = 



-0.880722 and = 50, E^c = -0.880976. 



It is worthwhile to note that an alternative summation for the first contribution in the 
right hand side of Eq.(8) has been derived in the appendix of ref. [IS] (cf. Eq.(A4)). In 
notations of the present paper, it becomes 

m=-oo m=i y (27rm)^ + L^K^ 

KLy 

(31) 

From this equation, we recover easily Eq.(lO) as k — > ; if k 7^ 0, the summation over m 
in the right hand side of Eq.(31) is slowly convergent and a lot of contributions have to 
be included to obtain a good accuracy, especially if z is small ; in agreement with Figures 
3 and 4. 



4 Conclusion. 

In view of the numerical tests performed in the previous section and since Lekner sums 
are not more efficient than Ewald sums : for Yukawa interaction potential in quasi-two 
dimensional systems, it is highly recommended to use the Ewald sums derived in ref.[T]. 
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List of Tables 



Table 1 : Numerical tests of the Coulomb limit of the Lekner sums for the Yukawa 
potentials. The configuration of the pair of particles is defined by {xi2, 2/12, -212) particle 
1 with a charge Qi = +1 located at (0,0,0) and particle 2 with a charge Q2 = — 1 
located at {xi2,yi2, zu)- The values of particle-particle Coulomb energies are extracted 
from original works by others. These energies have been obtained with the Lekner sums 
for Coulomb potential with truncations indexes {pc, fric) indicated in the Table and they 
agree with the values computed in ref. [5] with pc = 30 and rric = 3. Evaluations of 
particle-particle interaction energies in columns Direct, refer to evaluations performed by 
using straightforwardly the Lekner sums for Yukawa potentials given by Eqs.(l,8) and 
(9), and with pc = 50 and as indicated in the Table. The term Modified refers to 
evaluations of the energies of pairs with Eq.(28) that uses Taylor's expansions, and with 
Pc = 50 and rric as indicated in the Table (see also Figure 3 and 5). In the column Ewald, 
we give the values of the energy obtained for k = 1. by using Ewald sums with a = 6. 
and A; X A; = 16 X 16 (see ref. [1] for the numerical tests of the Ewald sums for Yukawa 
potentials). 
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Configurations 




Coulomb 




Direct 


Modified 


Direct 


Modified 


Ewald 


and References 




Limit 




K = 10-6 


K = 10-6 


K=l. 


K = 1. 


K= 1. 


{Xl2,yi2,zi2) 


















(0.5,0.5,0.) 


analytical 


-2.28472 


3 


-2.14344 


-2.28472 


-1.41986 


-1.29568 


-1.4307 


mm 






50 


-2.27482 




-1.43067 






(0.1,0.1,0.1) 


(14,3) 


-5.77211 


3 


-5.77212 


-5.77212 


-4.87007 


-4.77663 


-4.87041 








50 


-5.77212 




-4.87041 






(0.,0.,0.25) 


(10,3) 


-3.72483 


3 


-3.76028 


-3.72483 


-2.98415 


-2.9022 


-2.98362 


112J 






50 


-3.72731 




-2.98362 






(-0.25,-0.15,-0.2) 


(10,3) 


-2.82156 


3 


-2.83153 


-2.82157 


-2.04423 


-1.94915 


-2.04483 


I12j 






50 


-2.82226 




-2.04483 






(0.4,0.4,0.1) 


not 


-2.28608 


3 


-2.20082 


-2.28609 


-1.44881 


-1.32721 


-1.45555 


113J 


reported 




50 


-2.28015 




-1.45555 







Table 1 M. Mazars, Yukawa II : Lekner sums. 



18 



List of Figures 



Figure 1: Representation of AEi2{K,;Pc,mc) as contour plot for several values of k 
and pc- For the six figures, the particle 1 is located at (0,0,0) and particle 2 at {x,y,z) ; 
the common parameters for all figures are : Qi = Q2 = +1, = Ly = 20, —5 < x < 5, 
—5 < y < 5, 2; = 0.1 and rric = 3. For each figure, 20 isopotentials are represented and 
the values of AE12 are contained between -0.5 and 0.5 (i.e. isopotentials are separated by 
0.05 in reduced energy unit). For each figure, the value of k and Pc are indicated on the 
figures. For | x |> 5 or | y |> 5, computations have not been done because of irrelevant 
complications induced by the minimum image convention. 

Figure 2: Surface representation of AEi2{K;Pc,Tnc) with the same parameters as in 
Figure 1. The values of AE12 are contained between -0.5 and 0.5. 

Figure 3: Representation of ii^cc(Yukawa; n) for four configurations of a pair of Yukawa 
particles as functions of n. i?cc (Yukawa; k) are computed by using equations (8) and (9) 
inserted in (1), with pc = 50 and rric = 4. The particle 1 carries a charge Qi = +1 
and is located at (0,0,0) ; the particle 2 carries a charge Q2 = —1 and four positions of 
the particle 2 are considered. These configurations have already been considered in some 
previous works [5l [121 [13] (see also Table 1). For each configuration, the value of the 
Coulomb limit obtained from Eq.(l) with p^ = 50 and = 4 is indicated by a thick red 
horizontal line and limiting values are explicitly given. The differences with the values 
obtained by using Lekner sums for non-screened Coulomb interactions are given in Table 
1. 

Figure 4: Representation of i?cc(Yukawa; k) computed by using equations (8) and 
(9) inserted in (1) for a pair of Yukawa particles as function of for three values of 
Pc- The configuration of the pair is defined by : Particle 1, Qi = +1 located at (0,0,0) ; 
Particle 2, Q2 = —1 located at (0.4,0.4,0.1) [5l[13]. The value of the screening parameter 
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is K = 0.0001, close to the Coulomb limit. 

Figure 5: Representation of £'cc(Yukawa; k) as function of k, for the configuration 
(0.4,0.4,0.1) of particle 2. Curves with several values of rric in Eq.(l) and for a value 

rric = 3 in Eq.(28), that allows to overcome the lost of accuracy of Eq.(l) for small k, are 
represented ; for each curve Pc = 50. The inset shows the values obtained for 0.1 < k < 1 
where Taylor's expansions given in Eqs.(25) and (26) become inaccurate. 
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Figure 1: M. Mazars, Yukawa II : Lekner sums. 
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M. Mazars, Yukawa II 



Lekner sums. 
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Figure 3: M. Mazars, Yukawa II : Lekner sums. 
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Figure 4: M. Mazars, Yukawa II : Lekner sums. 
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Figure 5: M. Mazars, Yukawa II : Lekner sums. 
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